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Abstract 

Detailed  thermodynamic,  kinetic,  geometric,  and  cost  models  are  developed,  implemented,  and  validated  for  the  synthesis/design  and  operational 
analysis  of  hybrid  SOFC-gas  turbine-steam  turbine  systems  ranging  in  size  from  1.5  to  lOMWe.  The  fuel  cell  model  used  in  this  research  work 
is  based  on  a  tubular  Siemens-Westinghouse-type  SOFC,  which  is  integrated  with  a  gas  turbine  and  a  heat  recovery  steam  generator  (HRSG) 
integrated  in  turn  with  a  steam  turbine  cycle.  The  current  work  considers  the  possible  benefits  of  using  the  exhaust  gases  in  a  HRSG  in  order  to 
produce  steam  which  drives  a  steam  turbine  for  additional  power  output.  Four  different  steam  turbine  cycles  are  considered  in  this  research  work: 
a  single-pressure,  a  dual-pressure,  a  triple  pressure,  and  a  triple  pressure  with  reheat.  The  models  have  been  developed  to  function  both  at  design 
(full  load)  and  off-design  (partial  load)  conditions.  In  addition,  different  solid  oxide  fuel  cell  sizes  are  examined  to  assure  a  proper  selection  of 
SOFC  size  based  on  efficiency  or  cost.  The  thermoeconomic  analysis  includes  cost  functions  developed  specifically  for  the  different  system  and 
component  sizes  (capacities)  analyzed.  A  parametric  study  is  used  to  determine  the  most  viable  system/component  syntheses/designs  based  on 
maximizing  total  system  efficiency  or  minimizing  total  system  life  cycle  cost. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  combination  of  SOFCs  with  gas  turbines  is  one  of  the 
most  promising  power  generating  applications  [1,2].  The  SOFC 
stack  forms  the  combustor  unit  in  a  gas  turbine  system.  Com¬ 
pressed  air  is  fed  into  the  SOFC  stack  where  fuel  is  injected 
and  electrical  power  drawn  off.  Operating  close  to  a  46%  con¬ 
version  of  fuel  to  electrical  power  [3],  this  SOFC  then  provides 
pressurized  hot  gases  to  a  turbine  operating  at  35%  electrical  effi¬ 
ciency.  The  theoretical  overall  electrical  conversion  efficiency  of 
this  system  can  approach  65+%,  which  can  be  further  improved 
by  adding  a  steam  turbine  cycle  to  drive  the  overall  electrical 
efficiency  into  the  mid  seventies  [1,2]. 

The  objective  of  this  research  work  is  to  make  a  rigorous 
investigation  of  the  design  and  performance  characteristics  of 
hybrid  system  configurations  consisting  of  a  SOFC,  gas  tur¬ 
bine,  and  steam  turbine  for  stationary  power  applications  which 
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provide  power  to  a  large  number  of  residential/commercial 
buildings.  For  example,  a  lOMWe  hybrid  system  can  fulfill  the 
needs  of  2000  family  residences  based  on  an  average  four  person 
family  residence  in  the  US  which  requires  on  average  5  kWe. 
To  model  and  then  analyze  the  hybrid  system  configurations 
as  realistically  as  possible,  detailed  system/component  thermo¬ 
dynamic,  kinetic,  geometric,  and  cost  models  are  developed, 
implemented,  validated  and  then  used  to  conduct  a  parametric 
analysis  of  the  key  system/component  parameters  to  investigate 
both  thermodynamically  (efficiency  maximization)  and  eco¬ 
nomically  (total  life  cycle  cost  minimization)  the  advantages 
that  such  hybrid  systems  might  have  over  conventional  GT-ST 
combined  cycle  systems,  standalone  SOFC  systems,  and  hybrid 
SOFC-GT  systems. 

2.  System  layout 

As  a  starting  point  for  this  research  work,  the  modeling,  com¬ 
puter  code,  and  optimization  results  of  a  previously  developed 
1.5  MWe  hybrid  SOFC-GT  plant  Calise  et  al.  [4-6]  are  used. 
In  this  work,  the  gas  turbine  exhaust  mixture  is  re-circulated 
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Nomenclature 

A 

heat  transfer  area  (m2) 

D 

diameter  (m) 

L 

length  (m) 

m 

mass  flow  rate  (kg  s-1) 

^cell 

number  of  cells 

Q 

heat  transfer  rate  (kW) 

T 

temperature  (K) 

V 

volume  (m3) 

w 

power  output  (kW) 

Subscripts 

cell 

SOFC  stack 

comp 

compressor 

cond 

condenser 

cool 

cool  water  stream 

gas 

gas  side  stream 

GT 

gas  turbine 

HEC 

counter-flow  heat  exchanger 

i 

inlet  stream 

0 

outlet  stream 

PR 

pre-reformer 

pump 

pump 

ST 

steam  turbine 

steam 

steam  side  stream 

and  used  to  preheat  the  input  air  and  fuel  streams  by  means 
of  heat  exchangers,  while  the  remaining  energy  is  recovered 
to  heat  water  for  residential  usage,  while  in  the  current  work, 
this  system  has  been  modified  and  expanded  to  include  a  sec¬ 
ond  bottoming  cycle,  utilizing  various  types  (based  on  pressure 
level)  of  heat  recovery  steam  generators  and  a  steam  turbine.  In 
fact,  four  configurations  are  modeled  and  analyzed  in  detail  here 
with  the  variations  occurring  with  regard  to  the  steam  turbine 
(ST)  bottoming  cycle,  i.e.  a  single-pressure  level,  a  dual-pressure 
level,  and  a  triple-pressure  level  with  and  without  reheat. 
The  purpose  of  using  multiple-pressure  levels  is  to  achieve  a 
higher  power  output  from  the  steam  turbine  at  the  expense,  of 
course,  of  extra  equipment.  The  operation  of  the  SOFC-GT 
topping  cycle  (see  Fig.  1,  top  part)  can  be  summarized  as 
follows: 

•  Air  is  compressed  by  the  air  compressor  (AC)  up  to  the  fuel 
cell  operating  pressure.  The  air  is  then  brought  to  the  cathode 
inlet  of  the  SOFC  stack  (state  point  18).  Similarly,  fuel  is 
compressed  by  the  fuel  compressor  (FC)  and  then  brought  to 
the  anode  compartment  of  the  stack  (state  point  1). 

•  Both  fuel  and  air  can  by-pass  the  fuel  cell,  i.e.  a  certain  amount 
of  fuel  can  flow  directly  to  the  combustor  (C)  by-passing 
the  electrochemical  reaction  occurring  within  the  stack  (state 
point  23),  while  excess  air  can  flow  to  the  GT  (state  point  20). 

•  At  the  stack,  fuel  (state  point  24)  is  mixed  with  the  anode  re¬ 
circulation  stream  (state  point  5)  in  order  to  support  the  steam 
reforming  reaction  in  the  pre-reformer  and  in  the  anode  com¬ 


partment  of  the  fuel  cell.  The  mixture  at  state  point  25  consists 
of  methane  and  steam.  Thus,  in  the  pre-reformer  (PR),  the  first 
step  in  the  fuel  reforming  process  occurs.  The  energy  required 
to  support  the  pre-reforming  reaction  is  derived  from  the  hot 
stream  at  state  point  26.  The  non-reacted  fuel  at  state  point  2 
is  involved  in  the  internal  reforming  reaction  within  the  anode 
compartment  of  the  SOFC  stack.  Here,  it  is  converted  into  the 
hydrogen  that  participates  in  the  electrochemical  reaction. 

•  On  the  cathode  side,  air  is  first  preheated  by  a  counter-flow 
heat  exchanger  air  injection  pipe  (HEC)  and  then  brought 
into  the  annulus  (air  pipe)  of  the  SOFC  where,  at  the  three- 
phase  boundaries,  the  cathode  electrochemical  reaction  occur 
[1-2, 7, 8]. 

•  The  electrochemical  reactions,  occurring  in  the  fuel  cell,  pro¬ 
duce  DC  electrical  current  and  release  thermal  energy.  The 
first  of  these  is  converted  into  AC  current  by  the  inverter;  the 
latter  is  used  by  the  internal  reforming  reaction  and  to  heat  up 
the  fuel  cell  stack. 

•  The  high  energy  flow  rate  at  state  point  8  is  first  used  to  preheat 
air  in  the  counter-flow  heat  exchanger  and  then  to  supply 
energy  to  the  pre-reforming  reaction.  This  stream  at  state  point 
21  enters  the  gas  turbine. 

•  The  expansion  in  the  GT  supplies  mechanical  power  which 
in  turn  is  converted  into  electric  power. 

•  The  operation  of  the  steam  turbine  bottoming  cycle,  e.g.,  the 
triple  pressure  with  reheat  variation  (see  Fig.  1,  bottom  part), 
can  be  summarized  as  follows: 

•  The  GT  exhaust  stream  (state  point  33)  flows  to  the  heat  recov¬ 
ery  steam  generator  (HRSG).  The  gas  mixture  side  of  the 
HRSG  passes  through  the  ten  heat  exchanger  sections  -  high- 
pressure  (HP)  superheater  (SU),  reheater  (RH),  HP  evaporator 
(EV),  HP  economizer  (EC),  intermediate-pressure  (IP)  SU,  IP 
EV,  IP  EC,  low-pressure  (LP)  SU,  LP  EV,  and  LP  EC  -  and 
is  exhausted  at  state  point  34. 

•  The  superheated  steam  produced  by  the  HP  SU  (state  point 
35)  is  supplied  to  the  HP  stage  of  the  steam  turbine.  After 
expansion  the  cold  reheat  (state  point  64)  at  an  intermedi¬ 
ate  pressure  returns  to  the  HRSG  and  there  by  means  of  a 
reheater  is  superheated  (state  point  66)  and  returned  to  the 
IP/LP  steam  turbine  stage.  Also  the  IP  SU  (state  point  56) 
and  the  LP  SU  (state  point  48)  supply  superheated  steam 
to  the  double-admission  IP/LP  ST  which  during  expansion 
produces  mechanical  power  which  in  turn  is  converted  into 
electric  power  in  a  generator.  A  small  fraction  of  super¬ 
heated  steam  at  low  pressure  is  extracted  (state  point  37) 
to  the  deaerator  (DE)  to  be  used  later  on  for  feedwater 
preheating. 

•  The  wet  steam  (state  point  38)  is  then  condensed  in  the  con¬ 
denser  (CON).  The  condensate  (state  point  39)  enters  the 
condensate  pump  (CP)  and  is  then  pumped  to  the  DE  at  state 
point  40. 

•  In  the  DE,  any  air  oddments  and  impurities  contained  by  the 
water  are  removed  while  the  water  is  preheated  at  60  °C.  The 
preheated  water  (state  points  57,  49,  41)  enters  the  HP  FP 
(feedwater  pump),  IP  FP,  and  LP  FP,  and  is  then  pumped  to 
the  HP  EC,  IP  EC  and  LP  EC  at  state  points  58,  50,  and  42, 
respectively. 
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Fig.  1.  SOFC-GT  integrated  with  a  triple  pressure  with  reheat  ST  cycle. 


•  In  the  economizers,  water  is  heated  up  to  the  saturated  liquid 
point.  Then  it  is  evaporated  at  constant  temperature/pressure 
in  the  evaporators. 

•  Water  and  saturated  steam  are  separated  in  the  drums,  and  the 
steam  is  supplied  to  the  superheaters  where  it  is  superheated 
to  the  desired  live  steam  temperatures  and  fed  to  the  ST  to 
repeat  the  cycle. 

3.  Plant  model 

The  thermodynamic,  geometric,  kinetic  and  cost  models  were 
developed  in  MATFAB® .  They  are  used  to  simulate  the  behavior 
of  the  hybrid  fuel  cell  system  configuration.  The  model  is  based 
on  the  following  assumptions:  one-dimensional  flow;  steady 
state;  no  gas  leakage;  negligible  heat  losses  to  the  environment; 
negligible  kinetic  and  gravitational  terms  in  the  energy  balances. 

3.1.  Internal  reforming  SOFC  model 

The  overall  voltage  of  the  single  cell  can  be  calculated  as  a 
function  of  current  density,  temperatures,  pressures,  chemical 
composition,  and  geometric/material  characteristics  by  calcu¬ 
lating  the  difference  between  the  reversible  potential  and  all  the 
overvoltages  [6],  i.e. 

(1) 


where  V  is  the  actual  fuel  cell  potential,  E  the  open  circuit 
reversible  voltage  calculated  on  the  basis  of  the  Nernst  equa¬ 
tion,  yact  the  activation  overvoltage  calculated  on  the  basis  of 
the  Butler- Volmer  equation  and  experimental  correlations  for  the 
anode  and  cathode  exchange  current  densities,  V0hm  the  ohmic 
overvoltage,  and  Vconc  the  concentration  overvoltage.  Eq.  (1) 
suggests  that  in  the  case  of  SOFCs,  it  is  possible  to  neglect 
crossover,  fuel,  and  internal  current  losses.  The  assumptions 
and  calculation  details  for  all  aforementioned  overvoltages  are 
discussed  more  extensively  in  [6]. 

3.2.  P re -reformer  model 

One  of  the  main  advantages  of  using  high  temperature  fuel 
cells  is  the  possibility  of  feeding  the  SOFC  with  natural  gas 
directly,  since  the  reforming  process  can  be  supported  inside 
the  stack  [1,7-10].  In  practice,  however,  a  pre-reforming  pro¬ 
cess  is  usually  necessary.  The  pre-reformer  unit  consists  of  a 
number  of  tubes  located  inside  a  shell  and  filled  with  a  particu¬ 
lar  catalyst  [11,12].  The  reformate  gas  flows  inside  these  tubes. 
Hot  gases,  coming  from  the  combustor,  flow  inside  the  shell 
external  to  the  tubes,  supplying  the  thermal  energy  needed  to 
support  the  process,  since  the  energy  provided  by  the  exother¬ 
mic  water-gas  shift  reaction  is  not  sufficient  for  the  endothermic 
demethanization  of  the  reforming  process  [11,12].  The  pre- 
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reformer  heat-exchange  model  is  simulated  taking  into  account 
both  equilibrium  and  the  kinetics  as  described  in  [6]. 

3.3.  Gas  turbine  cycle  model 

The  SOFC-GT  subsystem  utilizes  an  air  compressor,  fuel 
compressor,  and  gas  turbine.  The  air  compressor  and  gas  tur¬ 
bine  are  connected  together  with  a  single  shaft.  The  shaft  is 
also  connected  to  an  electric  generator  converting  the  mechani¬ 
cal  power  to  electrical  power.  Mass  flow  rates  and  rotor  speeds 
are  corrected  on  the  basis  of  their  inlet  conditions  according  to 
[4-6,13,14]. 


3.4.  Counter-flow  heat  exchanger  air  injection  pipe, 
catalytic  combustor,  mixer,  inverter,  electric  generator 
models 

A  counter-flow  tube-in-tube  heat  exchanger  is  required  in 
order  to  simulate  the  heat  transfer  in  the  air  injection  pipe 
between  the  air  flowing  through  the  fuel  cell  air  tube  and  the 
stream  coming  from  inside  the  stack  [1-2,5-8,10,14-16].  The 
heat  exchange  is  simulated  on  the  basis  of  existing  models  in  [5] 
and  improved  to  include  the  effects  of  pressure  drops  and  to  take 
into  account  the  dependence  of  the  thermophysical  and  transport 
properties  on  temperature.  The  details  of  the  counter-flow  heat 
exchanger  heat  transfer  model  are  given  in  [4] . 

The  combustor  burns  any  non-reacted  fuel  coming  out  of 
the  fuel  cell  and,  therefore,  produces  thermal  energy  for  use 
elsewhere  in  the  system.  The  hybrid  plant  makes  use  of  three 
mixers.  These  are  necessary  for  the  operation  and  the  regulation 
of  the  plant.  The  exerted  electric  signal  needs  to  be  conditioned 
before  usage,  converted  to  AC  current,  and  filtered  from  possi¬ 
ble  oscillations.  This  is  done  by  a  DC- AC  inverter.  Similarly, 
the  mechanical  energy  produced  by  the  gas  turbine  must  be  con¬ 
verted  to  electric  power.  This  conversion  is  accomplished  by  an 
electric  generator. 

3.5.  Steam  turbine  cycle  model 

The  main  components  of  the  steam  turbine  cycle  include  a 
steam  turbine  with  an  electric  generator;  a  heat  recovery  steam 
generator  which  includes  economizer(s),  evaporator(s),  super¬ 
heaters),  and  a  reheater  (triple-pressure  reheat  cycle  only);  a 
condenser  which  is  dimensioned  according  to  the  turbine  exit 
pressure  and  mass  flow  rate  as  well  as  ambient  conditions;  a 
deaerator  heated  by  steam  extracted  from  the  steam  turbine;  a 
condensate  pump;  and  one  to  three  feedwater  pumps. 

The  HRSG  model  calculates  the  live  steam  mass  flow  rates 
and  also  the  exhaust  gas  conditions  at  the  HRSG  exit.  In  addi¬ 
tion,  it  sizes  the  different  types  of  heat  exchangers  included  in 
the  HRSG.  Depending  on  the  HRSG’s  number  of  pressure  lev¬ 
els,  the  corresponding  live  steam  mass  flow  rates  are  calculated. 
The  water/steam  conditions  at  the  inlet  and  exit  of  every  heat 
exchanger  are  defined  either  directly  by  the  desired  live  steam 
conditions  or  indirectly  through  conditions  on  the  saturation 
curve  [17].  An  important  parameter  defining  the  heating  surface 
and  performance  of  the  HRSG  is  the  pinch  point.  The  pinch- 


point  temperature  is  the  difference  between  the  evaporator’s 
outlet  temperature  on  the  water/steam  side  and  the  inlet  tem¬ 
perature  on  the  exhaust  gas  side.  The  lower  the  pinch  point,  the 
more  heating  surface  is  required  and  the  more  steam  is  generated 
[18]. 

The  desired  live  steam  temperatures  and  pressures  are  fixed. 
The  evaporator  drum  pressure  can  be  determined  based  on  a 
7-10%  loss  from  the  live  steam  pressure.  The  pinch  points  are 
also  selected  and  fixed.  The  energy  balances  on  the  gas  and  steam 
sides  are: 


Gsu+EV  —  ^GTexhcp(^SUin  ^EVout)  (2) 

GsU+EV  =  ^STin(^SUout  “  ^EVin)  (3) 


The  heat  transfer  rate  is  determined  from  Eq.  (2),  and  since  the 
two  heat  transfer  rates  on  the  left  hand  side  of  each  equation  are 
equal  to  each  other,  Eq.  (3)  is  solved  for  the  live  steam  mass  flow 
rate.  Using  simple  energy  balances,  identical  to  the  preceding 
ones,  all  temperatures  and  heat  transfer  rates  can  be  calculated 
for  all  the  heat  exchangers. 

For  the  geometric  models  of  the  heat  exchangers,  both  the 
LMTD  and  effectiveness-NTU  methods  are  used  depending  on 
the  exchanger.  The  geometric  models  are  needed  for  determining 
off-design  behavior.  All  the  heat  exchangers  are  shell-and-tube 
since  they  are  the  appropriate  type  for  compact  heat  recovery 
steam  generators  [19].  The  necessary  equations  for  shell-and- 
tube  heat  exchangers  are  obtained  from  [19].  A  detailed  analysis 
of  all  the  heat  exchangers’  geometric  models  can  be  found  in 
[20]. 

The  axial-flow  steam  turbine  can  be  single,  dual,  or  triple 
admission  depending  on  the  HRSG’s  pressure  level.  Further¬ 
more,  in  the  triple-pressure  reheat  cycle  configuration,  it  is 
divided  into  two  sections:  a  high  pressure  (HP)  section  and 
an  intermediate/low-pressure  (IP/LP)  section.  In  this  particular 
configuration,  the  HP  section  is  supplied  with  live  steam  from  the 
superheater.  After  expansion,  the  steam  returns  to  the  reheater 
in  the  HRSG  where  it  is  superheated  and  supplied  to  the  IP/LP 
section  for  further  expansion.  Following  expansion,  the  exhaust 
is  fed  to  the  condenser.  All  the  configurations  include  extraction 
outlets  for  deaerating/preheating.  The  mass  and  energy  balances 
for  the  triple-pressure  reheat  cycle  are 
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where  m? Ip.  ,  rnSL.  ,  mfJ  ,  and  m?p.  are  the  mass  flow 
rates  of  the  HP,  RH,  IP,  and  LP  superheated  steam  entering  the 
steam  turbine,  respectively,  and  the  mass  flow  rate  of  the 

steam  after  expansion  in  the  HP  section  of  the  steam  turbine. 
For  off-design  purposes  steam  turbine  maps  are  used  in  order  to 
capture  the  effects  of  geometry  on  turbine  performance.  To  gen¬ 
erate  these  maps  for  different  size  turbines,  data  is  taken  from 
[21]. 
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The  steam  turbine  cycle  includes  a  condensate  pump  and  one 
to  three  feedwater  pumps  depending  on  the  number  of  HRSG 
pressure  levels.  Since  the  thermodynamic  states  in  the  inlet  are 
known  and  the  outlet  thermodynamic  states  can  be  fixed  as 
desired,  what  is  left  is  a  calculation  of  the  pump  power  con¬ 
sumed.  The  corresponding  mass  and  energy  balances  are  given 
by 


.  pump 

m-  =  m 


pump 


Wpump  —  l/hin(^out  ^in)] 


pump 


(6) 

(7) 


steam  turbine  extraction,  md™  the  mass  flow  rate  of  the  deaer¬ 
ated/preheated  water  exiting  the  deaerator.  An  energy  balance 
on  this  component  yields 


•  dear  dea  •  dear  dea  •  dea i  dea  _  n 

Pin TYl\ n  in  O 


(12) 


where  hd™,  hd™,  and  hd^  are  the  enthalpies  for  the  correspond¬ 
ing  mass  flow  rates. 


4.  Cost  model 


where  mm  is  the  non-pressurized  mass  flow  rate  entering  the 
pump,  fh0 ut  the  pressurized  mass  flow  rate  exiting  the  pump, 
ffpump  the  pump  work  rate  consumption,  and  hp"mp  and  hp^p 
are  the  corresponding  enthalpies  for  the  mass  flow.  Again,  for 
off-design  purposes,  pump  maps  are  rescaled  and  modified  from 
actual  pump  maps  found  in  the  literature.  For  the  condensate 
pump,  a  map  from  [22]  for  a  centrifugal  type  pump  is  used 
while  for  the  feedwater  pumps,  a  map  for  a  displacement  type 
pump  from  [23]  is  employed. 

The  condenser,  which  is  a  shell-and-tube  heat  exchanger, 
receives  wet  steam  from  the  steam  turbine’s  exhaust  and  con¬ 
denses  it  to  a  saturated  liquid.  In  the  condensing  process,  the 
temperature  and  pressure  are  kept  constant.  For  the  purposes  of 
this  study,  they  have  been  fixed  at  3 1  °C  and  0.045  bar  as  in  [18]. 
The  working  side  mass  balance  is: 


^cond  _  ^cond 
rn  in  —  mout 


(8) 


where  mf°nd  is  the  mass  flow  rate  of  the  wet  steam  entering  the 
condenser  and  m^d  the  mass  flow  rate  of  the  saturated  liquid 
exiting.  The  heat  rejected  to  the  cooling  water  is  found  from  an 
energy  balance  on  the  condensing  steam,  i.e. 


Gcond  =  [<"d(Cd  -  hZf)]  (9) 

where  Gcond  is  the  rejected  heat  transfer  rate  and  /^°nd,  h™d  are 
the  enthalpies  for  the  corresponding  mass  flow  rates.  The  cooling 
water  mass  flow  rate  can  be  calculated  by  an  energy  balance  on 
the  cooling  water  entering  and  exiting  the  condenser: 


mcw  — 


Gcond 


(Tew  ,out 


a)c. 


pew 


(10) 


where  mcw  is  the  mass  flow  rate  of  the  cooling  water,  TCW;in 
and  rcw?0ut  are  the  inlet  and  outlet  cooling  water  temperatures, 
respectively,  and  Cp?cw  is  the  average  cooling  water  specific  heat. 
The  LMTD  method  is  applied  to  the  geometric  analysis  of  the 
condenser  analyzed  in  detail  in  [20]. 

The  deaerator  removes  dissolved  gases  and  impurities  from 
the  condensate  by  keeping  it  in  a  reservoir  at  the  state  of  a  sat¬ 
urated  liquid  absorbing  heat  extracted  from  the  steam  turbine 
at  a  pressure  slightly  higher  than  the  deaerator  pressure.  The 
corresponding  mass  and  energy  balances  are 


:  m 


•  dea 


(ID 


where  md ™  is  the  mass  flow  rate  of  the  saturated  liquid  com¬ 
ing  from  the  condensate  pump,  md™  the  mass  flow  rate  of  the 


For  the  thermoeconomic  analysis  of  the  plant,  appropriate 
cost  functions  must  be  formulated  to  include  the  purchase  cost 
for  every  component,  the  capital  cost  per  annum,  the  operating 
cost  per  annum,  and  the  total  cost  per  annum.  The  expressions 
for  all  the  component  purchase  costs  are  summarized  in  detail 
in  Table  1,  while  those  for  the  capital,  operating,  and  total  costs 
per  annum  are  summarized  in  Table  2. 

For  the  gas  turbine,  the  cost  function  proposed  by  [24]  is  used. 
For  the  centrifugal  compressors  (air  and  fuel  compressors),  the 
corresponding  costs  are  calculated  by  interpolating  data  from  the 
manufacturers  as  a  function  of  the  maximum  power  required  and 
using  information  provided  by  [25].  For  the  counter-flow  heat 
exchanger,  the  capital  cost  is  determined  on  the  basis  of  a  cost 
function  from  [26]  while  the  cost  of  the  SOFC  stack  is  esti¬ 
mated  with  reference  to  market  studies  in  which  the  expected 
cost  for  the  case  of  a  significant  increase  in  production  volume 
is  assumed.  A  detailed  work  performed  by  [9]  relates  the  SOFC 
purchase  cost  to  the  active  area  and  the  operating  temperature. 
Furthermore,  the  electric  energy  produced  by  the  SOFC  must  be 
filtered  by  an  inverter,  whose  cost  is  not  negligible  and  should, 
therefore,  be  taken  into  account  [9].  The  SOFC  system  also  con¬ 
sists  of  a  pre-reformer,  whose  cost  is  calculated  on  the  basis  of  its 
catalysts  volume  and  the  finned  exchange  area  [1 1-12,26]  which 
in  turn  is  related  to  the  number,  diameter,  and  length  of  tubes. 
Thus,  based  on  these  references  and  updating  the  functions  with 
literature  data,  the  pre-reformer  component  cost  function  is  for¬ 
mulated  by  [4].  The  total  cost  for  SOFC  auxiliary  devices  such 
as  the  combustor,  mixers,  and  by-pass  valves  are  calculated  as  a 
fixed  percentage  (10%)  of  the  stack  cost. 

For  the  steam  turbine  cycle,  all  cost  equations,  except  that  for 
the  steam  turbine,  are  based  on  [17]  and  have  been  appropriately 
adjusted  for  inflation  by  using  [27].  For  the  steam  turbine,  the 
cost  function,  which  is  based  on  the  steam  turbine  power  output, 
is  developed  based  on  [28].  For  the  HRSG  (which  includes  the 
drum  and  piping  costs),  the  total  cost  is  composed  of  the  cost  for 
the  various  heat  exchangers,  the  piping,  the  gas  conduit,  and  the 
pump.  It  is  based  on  a  function  used  by  [29] .  The  total  cost  of  the 
heat  exchangers  is  formed  by  the  sum  of  the  cost  for  the  various 
heat  exchange  units  (e.g.,  HP  superheater,  HP  evaporator,  HP 
economizer,  reheater,  etc.)  indicated  by  the  index  i.  Also  the 
LMTD  correction  factor,  Kt  is  based  on  the  logarithmic  mean 
temperature  difference,  ATim?/,  while  fPi,  /^steam,  and  fTigas 
are  cost  correction  factors.  The  cost  functions  for  the  piping 
and  the  gas  conduit  include  the  factors  fPi,  fTi  ste am,  and  fTigiXS . 
The  factors  introduce  a  sensitivity  of  cost  to  pressure  as  well 
as  to  steam  and  exhaust  gas  temperature.  The  pressure  factor  is 
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Table  1 

Component  cost  models 


Variable  description 

Model  equation 

Cgt 

Gas  turbine  component  cost  ($) 

CGt  =  (-98.328  ln(tVGT)  +  1318.5)IVgt 

Ccomp 

Compressor  component  cost  ($) 

Ccomp  =  91562((Vcomp/ 445)0'67 

Chec 

Counter-flow  heat  exchanger  component  cost  ($) 

Chec  =  130(Ahec/0.093)0,78 

CsOFC 

SOFC  stack  component  cost  ($) 

CsOFC  —  (^cells^f^cellCccll)(2.96Tcen  —  1907) 

Cinv 

Inverter  component  cost  ($) 

Cinv  =  105(WW500)ojo 

CpR 

Pre-reformer  component  cost  ($) 

Cpr  =  130(ApR,fln/0.093)0-78  +  3240(VPr)0-4  +  21,280.5VPR 

CaUx,SOFC 

SOFC  auxiliary  components  cost  ($) 

CaUx,SOFC  =  O.lOCsOFC 

Cst 

Steam  turbine  component  cost  ($) 

CSt  =  3644.3(WSt)°'7  -  61.3(WST)0'95 

fpi 

Heat  exchanger  pressure  factor 

fPi  =  0.097  Hpi /30)  +  0.9029 

/t, steam 

Steam-side  temperature  factor 

f%  steam  =  1  +  exp(Toutsteam  —  830/500) 

/T,gas 

Gas-side  temperature  factor 

/fgas  =  1  +  exp(rout,gas  —  990/500) 

Ki 

LMTD  correction  factor  (kWK-1) 

Ki  =  (Qi/ATlmJ) 

ChE(HRSG) 

HRSG’s  heat  exchangers  component  cost  ($) 

Che(hrsg)  =  365°y \fPi  f 7^steam  /t- gas  K0S). 

fpj 

Piping  pressure  factor 

fp.  =  0.097  l(/7y /30)  +  0.9029 

Cpiping 

HRSG’s  piping  component  cost  ($) 

Cpiping  —11?  820^  steam) 

Cgas 

HRSG’s  gas  conduit  cost  ($) 

j 

Cgas  =  658mJ£ 

Chrsg 

HRSG  component  cost  ($) 

Chrsg  =  Che(hrsg)  +  Cpiping  +  Cgas 

Ccond 

Condenser  component  cost  ($) 

Ccond  —  248Acond  +  659mcooi 

fr] 

Efficiency  correction  factor 

Fjj  =  1  +  (1  —  0.8/1  —  T^pump) 

Cpump 

Pump  component  cost  ($) 

Cpump  =  442(Wpump)°'711.41/, 

calculated  as  a  function  of  live  steam  pressure  P[  and  comes  from 
curve  fit  data  for  heat  exchangers  found  in  [26] .  The  temperature 
factors  are  developed  using  the  [29]  form  of  the  temperature 
correction  factors  and  the  fact  that  the  investment  of  superheaters 
is  about  twice  as  high  as  the  investment  cost  for  evaporators  [  17] . 
The  temperature  values  indicating  technical  limits  are  taken  from 
[17]. 

For  the  condenser,  the  cost  function  is  based  on  [29].  It  is 
calculated  as  a  function  of  the  condenser  surface  area,  Aconc i, 
and  the  cooling  water  mass  flow  rate,  mcoo l-  For  the  deaerator, 
the  cost  function  is  formulated  using  a  cost  function  found  in 
[26].  The  cost  function  for  the  pumps  is  taken  from  [29]  and 
calculates  the  cost  as  a  function  of  the  electric  power  consumed, 
Wpump,  and  an  efficiency  correction  factor/^. 

The  capital  cost  must  be  placed  on  an  annual  basis  in  order 
to  account  for  the  cost  of  the  investment  required.  This  annual 
cost  is  composed  of  the  depreciation  cost,  Cdep,  interest  on  the 
investment,  Cint,  maintenance  cost,  Cma i,  insurance  cost,  Cins, 
and  tax  cost,  Ctax-  The  depreciation  cost  is  based  on  the  fact  that 
the  equipment  deteriorates  with  time  [30]  and,  thus,  looses  value. 


This  loss  of  value  needs  to  be  distributed  over  the  lifetime  of  the 
component.  This  results  in  a  realistic  estimation  of  the  cost  of  the 
equipment  and  indicates  how  much  money  has  to  be  spent  every 
year  in  order  to  save  money  for  future  replacement  or  to  pay 
back  loans  if  the  equipment  was  purchased  with  outside  capital. 
In  the  context  of  thermoeconomic  modeling,  the  common  linear 
depreciation  method  is  used  for  this  cost  estimation.  Therefore, 
the  annual  depreciation  cost  is  determined  by  dividing  the  total 
purchase  cost,  Cpur,  by  the  depreciation  time,  ^dep,  measured  in 
years.  For  this  research  work,  n^v  has  been  assumed  to  be  10 
years  [4]. 

The  purchase  or  capital  cost  must  also  be  financed  from 
outside  sources  such  as  bank  loans.  The  associated  interest  is 
considered  a  cost  [30].  For  the  current  research  work,  some 
simplifying  assumptions  are  made:  (i)  a  single  interest  rate  is 
assumed  for  the  cost  of  borrowed  capital  as  well  as  for  the 
opportunity  cost  of  having  invested  ones  own  capital  and  (ii) 
the  capital  cost  is  distributed  over  the  lifetime  or  depreciation 
time  of  the  plant.  The  interest  rate,  i,  is  assumed  to  be  0.0926 
[30]. 


Table  2 

Capital,  operating,  and  total  cost  models 


Variable  description 

Model  equation 

Cdep 

Depreciation  cost  ($  year-1) 

Cdep  —  Cpur/n^ep?  Cint  —  (CpUr/ndep)G 

Cint 

Interest  on  outside  capital  cost  ($  year-1) 

Cmai  —  (Cpur/  ^dep)/mai 

Cmai 

Maintenance  cost  (Syear-1) 

Cins 

Insurance  cost  ($  year-1) 

Cjns  —  (  Cpur /^dep) /ins? 

Ctax 

Taxation  cost  (Syear-1) 

Ctax  —  (Cpur/ ^depX/tax  _ 

CCap 

Capital  cost  ($  year-1) 

Ccap  —  C^p  T  Cint  T  Cmai  ~\~  Cins  T  Ctax?  CGpe  =  CfVf/Vh 

Cope 

Operating  cost  (Syear-1) 

Ctotal  —  Ccap  +  C0pe 

Ctotal 

Total  cost  (Syear-1) 
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The  maintenance  cost  may  vary  over  the  lifetime  of  an  instal¬ 
lation  as  the  equipment  degrades  and  depends  largely  on  the 
number  of  operating  hours,  the  frequency  of  shutdowns  and  star¬ 
tups,  and  the  operating  environment.  A  total  maintenance  cost 
for  the  above  suggests  annual  maintenance  expenses  on  the  order 
of  6%  of  the  annual  depreciation  cost  [11].  Thus,/mai  which  is 
the  maintenance  cost  factor  is  0.06.  Similarly,  the  insurance  and 
taxation  cost  factors  are  chosen  as  0.2  and  0.54%,  respectively 

[17] . 

The  operating  cost  per  annum  is  based  on  [4].  In  this  cost 
function,  Cf  is  the  cost  of  fuel  in  $Nm-3,  Vf  is  the  volumetric 
flow  rate  of  the  fuel  in  Nm3  h-1,  and  Ah  are  the  annual  hours  of 
operation.  The  latter  is  assumed  to  be  8760  h  year-1  [4].  Once 
both  the  annual  operating  cost  and  the  capital  cost  per  annum  are 
known,  the  total  cost  per  annum  becomes  the  sum  of  the  annual 
capital  and  operating  costs. 

5.  Model  validation 

The  SOFC-GT-ST  hybrid  system  needs  to  be  validated  in 
order  to  have  confidence  in  the  model  predictions.  The  vali¬ 
dation  procedure  helps  determine  the  degree  of  accuracy  of 
the  model  and  any  possible  mismatches  and  discrepancies. 
Each  subsystem  of  the  hybrid  system  is  validated  separately 
using  manufacturer’s  data  for  the  SOFC  model  [4-6]  and 
measured  data  from  the  literature  for  the  steam  turbine  cycle 
models. 

As  mentioned  previously,  the  steam  turbine  cycle  subsys¬ 
tem  models  were  modeled  based  on  [18].  Fortunately,  the  same 
source  includes  measured  data  for  the  performance  of  all  four 
steam  turbine  cycle  configurations.  This  data  includes  temper¬ 
atures,  pressures,  mass  flow  rates,  and  steam  turbine  power 
output.  Detailed  comparisons  of  these  data  can  be  found  in 
[20].  The  comparison  indicates  good  agreement  between  the 
measured  values  from  the  literature  and  those  calculated  from 
the  models.  The  few  configurational  differences  which  do  exist 
are  justified  by  the  differences  which  exist  between  the  cycles 
found  in  the  literature  and  the  model  configurations  used  here. 
These  differences  include  the  following:  simplification  of  the 
HRSG’s  heat  exchangers  from  dual  flow  to  single  flow,  a  differ¬ 
ent  chemical  composition  for  the  GT  exhaust  gases,  a  pressure 
loss  (2%)  for  every  heat  exchanger  in  the  HRSG  for  all  model 
calculations  different  from  that  in  [  1 8] ,  and  a  model  pump  power 
expenditure  calculation  somewhat  different  than  that  found  in 

[18] .  The  slightly  lower  live  steam  production,  and  in  effect  the 
slightly  lower  ST  power  output  for  all  model  configurations  as 
compared  with  [18],  is  mainly  caused  by  the  absence  of  a  GT 
cooler. 

A  second  validation  was  also  made  to  check  in  more  detail 
the  steam  production  using  results  more  suitable  to  fuel  cell 
hybrid  systems.  These  results  are  from  a  simple,  single-pressure 
HRSG.  These  were  found  in  [2].  The  model  and  literature  results 
comparison  can  be  found  in  [20] .  The  minor  differences  between 
the  literature  and  the  model  results  are  caused  by  the  assumptions 
of  a  constant  specific  heat  throughout  the  HRSG  in  the  literature 
results  as  well  as  of  no  heat  exchanger  pressure  losses  in  the 
literature  results. 


6.  Parametric  study:  results  and  discussion 

The  purpose  of  the  parametric  study  is  to  either  maximize 
the  average  efficiency  or  minimize  the  total  cost  of  the  hybrid 
SOFC-GT-ST  power  plant.  To  achieve  this,  a  systematic  vari¬ 
ation  in  the  values  of  a  number  of  key  decision  variables  as 
well  as  the  relative  sizes  of  the  SOFC,  GT,  and  ST  for  a  given 
hybrid  plant  size  must  be  made.  The  key  decision  variable  ranges 
and  initial  values  are  chosen  on  the  basis  of  typical  hybrid  fuel 
cell  systems  found  in  [4-6] .  The  fuel  utilization  factor  is  varied 
from  0.75  to  0.90  in  0.05  increments.  Values  below  0.75  are  not 
applicable  because  such  values  cause  an  increase  in  temperature 
beyond  the  maximum  possible  turbine  inlet  temperature  (TIT) 
and,  therefore,  a  coupling  point  of  the  air  compressor  and  the 
gas  turbine  cannot  be  reached.  In  addition,  lower  fuel  utilization 
factors  result  in  lower  efficiencies  since  fuel  is  converted  mainly 
in  the  combustor. 

The  steam-to-carbon  ratio  is  varied  from  2  to  3.5  in  0.5 
increments.  Values  below  2  are  not  included  in  order  to  avoid 
problems  of  carbon  deposition  on  the  anode  of  the  SOFC  stack 
as  reported  in  [1-2,8]. 

The  SOFC  operating  temperature  is  varied  from  950  to 
1100°Cin50°C  increments.  A  value  beyond  1 100  is  not  used 
because  it  would  exceed  the  operating  limit  of  the  SOFC.  Also, 
a  value  lower  than  950  °C  at  full  load  (i.e.  at  the  design  point)  is 
infeasible  since  the  minimum  part  load  (25%)  SOFC  operating 
temperature  is  not  high  enough  to  heat  the  SOFC  stack. 

The  SOFC  operating  pressure  is  varied  from  7  to  lObar  in 
1  bar  increments.  A  value  beyond  10  bar  is  not  used  because  the 
differential  pressure  between  the  anode  and  the  cathode  com¬ 
partments  of  the  SOFC  would  exceed  the  operating  limits  set  by 
the  simulation  process  as  described  in  [6]. 

Finally,  the  unit  cost  of  fuel  is  varied  as  follows:  0.1,  0.3,  0.6, 
0.9,  and  1.2.  The  initial  value  of  0.3  is  the  one  published  in  [31] 
for  the  year  2005. 

Three  different  hybrid  plant  sizes  are  considered:  1.5,  5,  and 
lOMwe,  where  two  different  SOFC  sizes  are  considered  based 
on  current  density.  The  operating  current  density  for  the  selected 
fuel  cell  operates  from  100  to  650  mA  cm-2  [32].  Since  the  cur¬ 
rent  density  decreases  at  off-design,  the  small  SOFC  is  selected 
based  on  a  maximum  possible  current  density  of  550  mA  cm-2 
(full  load  or  design  condition)  while  the  large  SOFC  is  selected 
based  on  a  minimum  possible  current  density  of  100  mA  cm-2 
(25%  full  load  condition).  The  larger  SOFC  yields  higher  effi¬ 
ciencies  as  compared  with  the  smaller  SOFC  but  the  latter  has  a 
lower  capital  cost  which  is  significant  since  the  SOFC  purchase 
cost  is  the  most  dominant  of  all  the  equipment  purchase  costs. 
In  the  parametric  study  conducted  here,  the  smaller  SOFC  mini¬ 
mizes  the  total  cost  while  the  larger  one  maximizes  the  efficiency. 
Also,  the  gas  turbine  and  steam  turbine  sizing  is  determined 
based  on  the  extra  power  needed  to  reach  the  desired  total  power 
output. 

The  performance  of  each  individual  system  is  analyzed  at 
full  and  part  load  conditions  to  determine  the  average  and  total 
efficiencies  and  total  operating  cost.  The  load  profile  is  based 
on  a  2-day  (one  average  winter  day  and  one  average  summer 
day),  electrical  power  demand  profile  for  an  average  four  person 
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family  household  [33]  scaled  appropriately  to  coincide  with  the 
three  different  sized  hybrid  plants  analyzed  here  and  extended 
over  an  entire  year.  The  average  efficiency  for  the  plant  becomes 
a  time-averaged  value  based  on  the  time  intervals.  It  is  defined 
as  follows: 


Table  3 


Time  intervals  for  Eq.  (13) 

Variable  description 

Value  (days) 

Variable  description 

Value  (days) 

A£i  25%  load 

137 

A f3  75%  load 

106 

A?2  50%  load 

46 

A  A  100%  load 

76 

7?  a ve  — 


A*i77i  +  At2T]2  +  +  A  £4.774 

EM 


(13) 


where  the  At\,  A £2,  A £3,  and  A £4  are  the  time  intervals  cor¬ 
responding  to  25,  50,  75,  and  100%  loads  (values  shown  in 
Table  3),  and  the  771 , 772, 773,  and  774  are  the  hybrid  plant  energetic 
electrical  efficiencies  corresponding  to  the  above  time  intervals. 


In  addition,  the  total  efficiency  is  defined  as 


__  EWiAti 

mot  E  Gh  A  ti 


(14) 


where  VV,  is  the  power  output  corresponding  to  25,  50,  75, 
and  100%  load,  and  <2h  the  heat  input  to  the  system  based 


Table  4 

Parametric  results  for  the  lOMWe  triple-pressure  w/RH  ST  hybrid  system 

Decision  variables  Objective  functions  Power  output  (MW) 


uf 

Tsofc 

S/C 

PSOFC 

cfuel 

7ave 

Qotal 

*7tot 

WSOFC 

Wqt 

WSX 

Hybrid  system  employed  with  a  small  SOFC 

0.75  1000  2  8 

0.3 

0.5864 

3,846,790 

0.6240 

7.863 

2.308 

1.337 

0.8 

1000 

2 

8 

0.3 

0.5852 

3,852,258 

0.6213 

7.902 

2.356 

1.398 

0.85 

1000 

2 

8 

0.3 

0.5889 

3,637,267 

0.6242 

7.958 

2.400 

1.454 

0.9 

1000 

2 

8 

0.3 

0.5685 

3,741,024 

0.6039 

8.001 

2.447 

1.510 

0.85 

950 

2 

8 

0.3 

0.5665 

3,690,959 

0.6030 

7.999 

2.439 

1.407 

0.85 

1000 

2 

8 

0.3 

0.5889 

3,637,267 

0.6242 

7.958 

2.400 

1.454 

0.85 

1050 

2 

8 

0.3 

0.5960 

3,660,018 

0.6300 

7.902 

2.358 

1.506 

0.85 

1100 

2 

8 

0.3 

0.6032 

3,692,544 

0.6368 

7.856 

2.310 

1.548 

0.85 

1000 

2 

8 

0.3 

0.5889 

3,637,267 

0.6242 

7.958 

2.400 

1.454 

0.85 

1000 

2.5 

8 

0.3 

0.5822 

3,653,267 

0.6178 

8.006 

2.447 

1.405 

0.85 

1000 

3 

8 

0.3 

0.5810 

3,644,914 

0.6166 

8.059 

2.489 

1.360 

0.85 

1000 

3.5 

8 

0.3 

0.5716 

3,689,910 

0.6079 

8.102 

2.537 

1.307 

0.85 

1000 

2 

7 

0.3 

0.5848 

3,645,623 

0.6192 

7.906 

2.354 

1.403 

0.85 

1000 

2 

8 

0.3 

0.5889 

3,637,267 

0.6242 

7.958 

2.400 

1.454 

0.85 

1000 

2 

9 

0.3 

0.5869 

3,641,883 

0.6214 

8.001 

2.441 

1.498 

0.85 

1000 

2 

10 

0.3 

0.5852 

3,648,621 

0.6197 

8.054 

2.496 

1.543 

0.85 

1000 

2 

8 

0.1 

0.5889 

1,744,829 

0.6242 

7.958 

2.400 

1.454 

0.85 

1000 

2 

8 

0.3 

0.5889 

3,637,267 

0.6242 

7.958 

2.400 

1.454 

0.85 

1000 

2 

8 

0.6 

0.5889 

6,415,426 

0.6242 

7.958 

2.400 

1.454 

0.85 

1000 

2 

8 

0.9 

0.5889 

9,217,784 

0.6242 

7.958 

2.400 

1.454 

0.85 

1000 

2 

8 

1.2 

0.5889 

12,032,142 

0.6242 

7.958 

2.400 

1.454 

Hybrid  system  employed  with  a  large  SOFC 

0.75  1000  2  8 

0.3 

0.6298 

3,798,306 

0.6616 

8.303 

2.154 

1.187 

0.8 

1000 

2 

8 

0.3 

0.6367 

3,756,709 

0.6659 

8.347 

2.203 

1.126 

0.85 

1000 

2 

8 

0.3 

0.6529 

3,697,652 

0.6836 

8.391 

2.254 

1.079 

0.9 

1000 

2 

8 

0.3 

0.6304 

3,781,486 

0.6605 

8.442 

2.304 

1.005 

0.85 

950 

2 

8 

0.3 

0.6447 

3,625,057 

0.6741 

8.440 

2.997 

1.007 

0.85 

1000 

2 

8 

0.3 

0.6529 

3,697,652 

0.6836 

8.391 

2.254 

1.079 

0.85 

1050 

2 

8 

0.3 

0.6412 

3,834,151 

0.6707 

8.344 

2.201 

1.130 

0.85 

1100 

2 

8 

0.3 

0.6386 

3,926,399 

0.6677 

8.301 

2.149 

1.189 

0.85 

1000 

2 

8 

0.3 

0.6529 

3,697,652 

0.6836 

8.391 

2.254 

1.079 

0.85 

1000 

2.5 

8 

0.3 

0.6433 

3,737,498 

0.6722 

8.442 

2.305 

1.011 

0.85 

1000 

3 

8 

0.3 

0.6357 

3,748,615 

0.6657 

8.496 

2.347 

0.952 

0.85 

1000 

3.5 

8 

0.3 

0.6337 

3,774,777 

0.6641 

8.534 

2.396 

0.894 

0.85 

1000 

2 

7 

0.3 

0.6368 

3,823,306 

0.6658 

8.346 

2.203 

1.008 

0.85 

1000 

2 

8 

0.3 

0.6529 

3,697,652 

0.6836 

8.391 

2.254 

1.079 

0.85 

1000 

2 

9 

0.3 

0.6348 

3,840,558 

0.6644 

8.435 

2.301 

1.136 

0.85 

1000 

2 

10 

0.3 

0.6316 

3,871,101 

0.6610 

8.493 

2.348 

1.189 

0.85 

1000 

2 

8 

0.1 

0.6529 

1,985,021 

0.6836 

8.391 

2.254 

1.079 

0.85 

1000 

2 

8 

0.3 

0.6529 

3,697,652 

0.6836 

8.391 

2.254 

1.079 

0.85 

1000 

2 

8 

0.6 

0.6529 

6,281,601 

0.6836 

8.391 

2.254 

1.079 

0.85 

1000 

2 

8 

0.9 

0.6529 

8,859,549 

0.6836 

8.391 

2.254 

1.079 

0.85 

1000 

2 

8 

1.2 

0.6529 

11,449,497 

0.6836 

8.391 

2.254 

1.079 
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on  the  LHV  and  fuel  flow  rates  for  the  aforementioned  time 
intervals. 

Results  for  all  the  hybrid  systems  modeled  and  simulated 
here  are  analyzed  thoroughly  in  [20].  A  representative  set  of 
results  for  the  parametric  study  performed  for  the  lOMWe 
SOFC-GT-ST  hybrid  system  with  a  triple  pressure  with  reheat 
ST  cycle  is  shown  in  Table  4  and  Figs.  2  and  3.  The  top  half  of 
Table  4  shows  the  results  for  the  small  SOFC,  while  the  bottom 
half  shows  those  for  the  large  SOFC.  The  objective  function  is 
the  average  efficiency  when  using  a  large  SOFC  and  the  total 
cycle  cost  for  the  small  SOFC.  Figs.  2  and  3  show  the  deci¬ 
sion  variable  values  plotted  against  the  corresponding  objective 
function  values. 


From  an  observation  of  the  trend  lines  shown  in  these  figures, 
some  general  remarks  can  be  made  that  apply  to  all  the  hybrid 
plants  analyzed  here.  The  optimum  fuel  utilization  factor  (Uf)  is 
0.85.  For  the  small  SOFC  where  the  cost  is  the  optimizing  func¬ 
tion  of  interest,  the  minimum  total  cost  is  achieved  at  this  value 
because  although  the  capital  cost  decreases  slightly  at  lower  val¬ 
ues  of  Uf,  the  more  dominating  operating  cost  is  decreased  at 
higher  values  of  Uf.  Therefore,  the  higher  efficiency  is  the  main 
reason  that  the  higher  fuel  utilization  is  more  economical  even 
though  only  slightly  so. 

For  the  large  SOFC,  where  the  average  efficiency  is  the  opti¬ 
mizing  function  of  interest,  as  expected,  the  higher  the  fuel 
utilization  factor  is  the  higher  the  efficiency.  The  optimum  value 
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Fig.  2.  Optimizing  variables  vs.  objective  functions  for  the  lOMWe  triple-pressure  w/RH  ST  small  SOFC. 
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Fig.  3.  Optimizing  variables  vs.  objective  functions  for  the  lOMWe  triple-pressure  w/RH  ST  large  SOFC. 


is  0.85  and  not  0.90  because  at  the  latter  value  although  the 
SOFC  efficiency  is  slightly  higher,  the  total  plant  efficiency 
drops  because  not  much  heat  is  left  for  recovery  by  the  gas 
turbine  and  the  steam  turbine  since  almost  all  the  hydrogen  pro¬ 


duced  by  the  internal  reforming  reactions  is  consumed  within 
the  fuel  cell  by  the  anode  electrochemical  reaction  and  more 
efficient  stacks  release  less  heat.  Therefore,  the  electrical  power 
produced  by  the  SOFC  increases,  causing  a  raise  in  its  electro- 
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chemical  rate  of  reaction.  On  the  other  hand,  this  effect  decreases 
significantly  the  turbomachinery  efficiencies. 

The  optimum  SOFC  operating  temperature  for  both  SOFC 
sizes  and  all  models  is  1000  °C.  For  the  small  SOFC,  where 
the  total  cost  is  the  optimizing  function,  the  lower  the  SOFC 
operating  temperature  the  lower  the  SOFC  capital  cost.  This 
trend  reaches  a  minimum  at  1000  °C  and  not  950  °C  because 
the  operating  cost  manages  to  slightly  overcome  the  capital 
cost.  For  the  large  SOFC,  where  the  average  efficiency  is  the 
optimizing  function,  the  best  SOFC  operating  temperature  is 
still  1000  °C  because  although  the  SOFC  stacks  operate  slightly 
more  efficiently  at  increased  temperatures,  the  air  and  fuel 
compressors  require  significantly  higher  power  consumption, 
thereby  decreasing  slightly  the  overall  efficiency. 

The  steam-to-carbon  ratio  reaches  its  optimum  efficiency  and 
lower  cost  at  the  lowest  possible  value  of  2.  The  capital  cost 
remains  constant  throughout  the  S/C  variation;  and,  therefore, 
the  only  significantly  varying  cost  is  the  operating  cost  which 
will  be  determined  by  the  overall  efficiency.  The  average  effi¬ 
ciency  slightly  decreases  as  the  S/C  increases  since  higher  steam 
partial  pressures  cause  higher  Nernst  overvoltages  [6].  Thus, 
higher  efficiencies  and  lower  operating/total  costs  are  achieved 
at  lower  S/C  values. 

Theoretically  increasing  the  SOFC  operating  pressure  will 
result  in  an  increase  of  the  cell  voltage  because  of  the  higher  reac¬ 
tant  partial  pressure  available,  therefore,  improving  the  SOFC 
efficiency.  On  the  other  hand,  the  expansion  of  the  gas  in  the 
turbomachinery  is  in  a  temperature  region  where  the  turbo¬ 
machinery  produces  less  power  as  the  pressure  increases  as 
indicated  in  [34] .  Also,  the  additional  thermal  input  to  the  sys¬ 
tem  for  delivering  air  and  fuel  at  the  desired  conditions  leads 
to  a  decreasing  trend  for  the  overall  efficiency  when  the  pres- 


Table  5 

Optimal  component  costs  for  the  lOMWe  triple-pressure  w/RH  ST  hybrid 
system 


Component 

Cost 

Large  SOFC 

Small  SOFC 

SOFC 

5,031,872 

2,334,223 

GT 

1,595,840 

1,662,272 

Combustor/mixer 

1,236,623 

941,679 

HRSG 

584,085 

660,761 

ST 

438,749 

536,340 

Inlet  air  tubes 

345,616 

189,093 

AC 

258,163 

264,960 

Deaerator 

57,278 

67,761 

Pre-reformer 

45,124 

36,994 

FC 

41,321 

42,856 

Condenser 

8,262 

8,262 

Pumps 

7,958 

9,734 

Total  investment 

9,650,891 

6,754,935 

sure  is  increased  beyond  a  certain  level.  Therefore,  the  optimum 
efficiency  for  both  the  total  cost  and  efficiency  objectives  is  at 
an  intermediate  point  (8  bar)  where  the  capital  cost  is  not  as 
high  as  compared  to  a  capital  cost  for  an  even  higher  pressure 
where  the  capital  cost  of  the  turbomachinery  would  increase  to 
accommodate  the  higher  pressure  needed.  On  the  other  hand, 
a  lower  operating  pressure  causes  lower  efficiencies  and,  there¬ 
fore,  although  the  capital  cost  is  lower  the  operating  cost  is  higher 
resulting  in  an  increased  total  cost. 

Finally,  the  unit  cost  of  fuel  variation  helps  determine  how 
efficiently  the  overall  hybrid  plant  should  be  designed.  At  lower 
costs  of  fuel,  e.g.,  0.1  and  0.3  $  Nm-3,  the  capital  cost  competes 
more  evenly  with  the  operating  cost  to  determine  the  optimum 


Table  6 

Cost  and  efficiency  breakdown  for  all  optimal  models 


Configuration 

Description 

Wsofc  (MW) 

Total  cost 

Operating  cost 

Capital  cost 

7ave 

ritot 

7?  max 

1 

1.5  MWe  single-pressure  ST  small  SOFC 

1.275 

604,743 

438,971 

165,772 

0.5696 

0.6022 

0.6638 

2 

1.5  MWe  single-pressure  ST  large  SOFC 

1.330 

613,655 

405,328 

208,327 

0.6235 

0.6525 

0.6971 

3 

1.5  MWe  dual-pressure  ST  small  SOFC 

1.275 

602,033 

435,146 

166,887 

0.5705 

0.6031 

0.6647 

4 

1.5  MWe  dual-pressure  ST  large  SOFC 

1.330 

611,026 

401,575 

209,451 

0.6243 

0.6533 

0.6978 

5 

1.5  MWe  triple-pressure  ST  small  SOFC 

1.238 

599,321 

431,428 

167,893 

0.5712 

0.6035 

0.6654 

6 

1.5  MWe  triple-pressure  ST  large  SOFC 

1.314 

608,046 

397,831 

210,215 

0.6251 

0.6536 

0.6987 

7 

1.5  MWe  triple  RH-pressure  ST  small  SOFC 

1.234 

598,691 

430,598 

168,093 

0.5716 

0.6038 

0.6657 

8 

1.5  MWe  triple  RH-pressure  ST  large  SOFC 

1.300 

607,446 

397,031 

210,415 

0.6255 

0.6539 

0.6992 

9 

5  MWe  single-pressure  ST  small  SOFC 

4.116 

1,887,264 

1,437,339 

449,925 

0.5791 

0.6131 

0.6795 

10 

5  MWe  single-pressure  ST  large  SOFC 

4.333 

1,935,053 

1,340,421 

594,632 

0.6264 

0.6573 

0.7115 

11 

5  MWe  dual-pressure  ST  small  SOFC 

4.069 

1,872,891 

1,411,966 

460,925 

0.5808 

0.6149 

0.6815 

12 

5  MWe  dual-pressure  ST  large  SOFC 

4.325 

1,920,649 

1,316,017 

604,632 

0.6279 

0.6589 

0.7132 

13 

5  MWe  triple-pressure  ST  small  SOFC 

4.054 

1,871,891 

1,405,966 

465,925 

0.5814 

0.6155 

0.6825 

14 

5  MWe  triple-pressure  ST  large  SOFC 

4.317 

1,924,649 

1,310,017 

614,632 

0.6286 

0.6596 

0.7142 

15 

5  MWe  triple  RH-pressure  ST  small  SOFC 

4.037 

1,869,891 

1,402,966 

466,925 

0.5817 

0.6159 

0.6837 

16 

5  MWe  triple  RH-pressure  ST  large  SOFC 

4.281 

1,923,649 

1,307,017 

616,632 

0.6289 

0.6601 

0.7157 

17 

10  MWe  single-pressure  ST  small  SOFC 

7.999 

3,630,071 

2,831,362 

798,710 

0.5876 

0.6222 

0.6922 

18 

10  MWe  single-pressure  ST  large  SOFC 

8.495 

3,717,621 

2,603,916 

1,113,705 

0.6473 

0.6766 

0.7272 

19 

10  MWe  dual-pressure  ST  small  SOFC 

7.929 

3,630,067 

2,819,358 

810,710 

0.5880 

0.6228 

0.6943 

20 

10  MWe  dual-pressure  ST  large  SOFC 

8.485 

3,717,624 

2,591,919 

1,125,705 

0.6478 

0.6772 

0.7292 

21 

10  MWe  triple-pressure  ST  small  SOFC 

7.905 

3,640,067 

2,805,358 

834,710 

0.5884 

0.6234 

0.6961 

22 

10  MWe  triple-pressure  ST  large  SOFC 

8.474 

3,715,624 

2,577,919 

1,137,705 

0.6480 

0.6774 

0.7292 

23 

10  MWe  triple  RH-pressure  ST  small  SOFC 

7.958 

3,606,144 

2,822,572 

783,573 

0.5891 

0.6242 

0.6973 

24 

10  MWe  triple  RH-pressure  ST  large  SOFC 

8.391 

3,697,652 

2,577,948 

1,119,703 

0.6529 

0.6836 

0.7370 
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Fig.  4.  Cost  breakdown  for  all  optimal  configurations. 


system.  On  the  other  hand,  at  higher  values  of  unit  cost  of  fuel, 
e.g.,  0.9  and  1.2  $Nm-3,  the  operating  cost  increases  far  and 
away  above  the  capital  cost  and,  therefore,  a  minimization  of 
fuel  consumption  is  required.  In  such  a  case,  a  more  efficient 
system  is  required,  and,  therefore,  the  large  SOFC,  although 
having  a  higher  capital  cost  than  the  smaller  SOFC  would  be 
selected. 

The  component  cost  breakdown  for  the  lOMWe  triple¬ 
pressure  reheat  ST  hybrid  system  is  shown  in  Table  5.  As 
expected,  the  most  dominant  component  cost  is  the  SOFC 
purchase  (depreciation)  cost  (especially  when  using  the  large 
SOFC),  while  the  turbomachinery,  combustor,  and  HRSG  costs 
are  also  significant. 


All  twenty-four  system  configurations  are  compared  in  terms 
of  cost  (total,  operating,  and  capital)  and  efficiency  (maximum, 
total,  and  average)  in  Table  6.  From  Fig.  4  it  can  be  concluded  that 
as  expected  the  higher  the  degree  of  complexity  in  the  heat  recov¬ 
ery  steam  generator,  the  higher  the  capital  cost  is.  On  the  other 
hand,  the  operating  cost  decreases  as  this  complexity  increases 
because  of  the  increasing  efficiency.  This  decrease  in  operating 
cost  is  significant  enough  to  balance  and  even  slightly  decrease 
the  total  cost  for  a  more  complex  system.  In  terms  of  SOFC 
size,  the  configurations  equipped  with  a  smaller  SOFC  also  as 
expected  have  a  total  cost  lower  than  those  with  the  larger  SOFC. 

Finally,  from  Fig.  5  it  is  evident  that  the  efficiency  is  higher 
for  a  system  equipped  with  a  larger  SOFC  than  those  with  a 
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Fig.  5.  Efficiency  breakdown  for  all  optimal  configurations. 
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smaller  SOFC  for  any  given  configuration  pair.  In  addition,  a 
global  comparison  of  all  twenty-four  systems  shows  that  the 
most  efficient  system  is  a  system  at  the  larger  power  capacity 
level,  i.e.  the  lOMWe  configuration  is  more  efficient  than  the 
5  MWe  configuration  or  the  5  MWe  configuration  is  more  effi¬ 
cient  than  the  1.5  MWe  configuration.  This  is  mainly  due  to  the 
higher  efficiencies  achieved  by  the  turbomachinery  (gas  turbine, 
steam  turbine,  air  compressor,  fuel  compressor,  pumps)  at  the 
larger  capacities. 

7.  Conclusions 

The  high  efficiencies  developed  by  some  of  the  hybrid  sys¬ 
tems  are  of  great  interest  since  they  show  the  potential  for 
exceeding  those  of  the  best  commercial  heat  engine  cycles 
currently  available  or  projected.  For  instance,  the  10  MWe 
SOFC-GT-ST  hybrid  triple  pressure  with  reheat  system  exhibits 
efficiencies  (maximum  efficiency  of  73.8%,  an  average  effi¬ 
ciency  of  65.3%,  and  a  total  efficiency  of  68.4%)  that  cannot 
be  matched  by  other  conventional  and  non-conventional  cycles 
(e.g.,  a  standalone  SOFC,  SOFC-GT  hybrid  cycles,  etc.).  Inter¬ 
estingly,  the  SOFC-GT-ST  system  develops  high  efficiencies 
at  off-design  conditions  as  well.  However,  the  off-design  strat¬ 
egy  followed  in  this  research  work  was  a  simplistic  one  which 
involved  the  lowering  of  the  fuel  flow  rate  while  keeping  con¬ 
stant  the  air  flow  rate.  This  strategy  was  necessary  because  a 
constant  air-to-fuel  ratio  strategy  (which  maintains  high  effi¬ 
ciencies)  has  a  very  restricted  field  of  operation  (above  80% 
of  full  load)  not  applicable  for  the  load  profile  (down  to  25% 
of  full  load)  considered  in  this  research  study.  Therefore,  this 
simplistic  strategy  creates  conditions  at  which,  for  example,  the 
SOFC  operating  temperature  drops  significantly,  leading  to  the 
difficulty  of  maintaining  higher  efficiencies  than  those  actually 
exhibited.  Thus,  a  more  in-depth  analysis  of  the  off-design  strat¬ 
egy  should  be  done  to  see  if  a  better  strategy  can  indeed  be  found 
(e.g.,  the  removal  of  part  of  the  SOFC  stack  at  lower  loads  could 
conceivably  maintain  higher  efficiencies). 

The  parametric  study  identified  a  number  of  unforeseen  com¬ 
plexities  which  only  became  evident  after  the  integration  and 
development  of  the  total  system  configurations.  These  difficul¬ 
ties  included  the  proper  selection  of  the  SOFC  stack  size  and 
the  difficulty  of  finding  the  proper  steam  turbines  to  match  the 
system.  For  a  realistic  system,  a  1.5  MWe  SOFC-GT-ST  is  not 
as  attractive  and  efficient  as  a  5  or  a  10  MWe  system  because  the 
gas  turbine  and  especially  the  steam  turbine  are  very  inefficient 
at  small  sizes  resulting  in  lower  overall  system  efficiencies. 

A  careful  selection  of  component  designs  in  the  steam  turbine 
cycle  was  made  to  achieve  efficient  conversion  of  the  thermal 
energy  to  power  output  based  on  the  thermodynamic,  geometric, 
and  cost  models.  The  uniqueness  of  the  system  required  in  many 
instances  the  rescaling  and  remodeling  of  existing  components 
to  fulfill  the  needs  of  the  system.  Thus,  if  not  careful,  one  risks 
making  unrealistic  selections  and  cost  analyses  of  equipment. 
Therefore,  the  design  of  a  hybrid  SOFC-GT-ST  power  plant 
must  focus  on  all  the  components  and  not  only  on  the  SOFC. 
Special  attention  must  be  given  in  the  coupling  of  the  turbo¬ 
machinery  with  the  SOFC  and  the  heat  exchangers  in  order  to 


achieve  the  maximum  benefits  offered  by  the  hybrid  system. 
The  air  compressor  has  an  unexpectedly  high  consumption  of 
power  (about  75%  of  the  gas  turbine  power  output)  meaning 
that  a  study  must  be  made  in  order  to  minimize  this  trend  if 
possible.  Based  on  current  technology,  the  SOFC  operates  at 
pressures  of  7-9  bars  meaning  that  a  high  degree  of  compres¬ 
sion  is  required  to  fulfill  this  need.  New  developments  in  lower 
temperature  and  pressure  SOFCs  may  benefit  the  overall  system 
since  less  compression  will  be  required. 

Since  the  SOFC-GT-ST  system  involves  a  large  amount  of 
equipment  with  a  much  larger  number  of  decision  variables 
than  actually  considered  in  this  parametric  study,  a  more  com¬ 
plete  optimization  of  the  systems  should  be  done  in  order  to 
determine  more  detailed  syntheses/designs  than  those  presented 
here. 

Finally,  since  SOFCs  are  not  fully  commercialized,  a  more 
accurate  economic  analysis  than  that  made  here  cannot  be  made 
at  this  time.  The  high  capital  cost  suggested  in  this  research 
work  (even  adjusted  for  production  volume)  could  decrease  in 
the  near  future  leading  to  minimized  cost  syntheses/designs 
which  exhibit  even  higher  efficiencies  than  those  determined 
here. 
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